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I will sketch a proof that the short distance behavior of the even scaling functions for the 

Q ■ Ising model that arise by taking the scaling limit of the Ising correlation functions from 

O , below the critical temperature is given by the Luther-Peschel formula [6] (see below). The 

! fact that the Luther-Peschel formula is consistent with conformal field theory insights 

^ I into the correlations for the large scale limit of the critical Ising model has lead to the 

conviction in the Physics community that this formula is right [4]. The critical scaling 

limit for the Ising model is obtained by sitting right at Tc (the critical temperature) and 

extracting the large scale limit of the correlations at this fixed temperature. However, as 

C/^ [ far as I am aware there is no mathematical proof that the large scale limit of the critical 

Ph I Ising model even exists (with the exception of the asymptotics of the diagonal two point 

correlation at the critical point [15]). In this note we will not address the question of the 

critical (massless) scaling limit for the Ising model but instead we will consider the massive 

scaling limit obtained by looking at the lattice correlations at the scale of the correlation 

>• I length as the temperature approaches the critical temperature. This scaling limit was first 

considered in [14] and [11] and a mathematical proof of its existence was established in [7]. 

O ' What we will sketch a proof of here is that the short distance asymptotics of this scaling 

^ . limit for the even scaling function from below Tc is given by the Luther-Peschel formula. 

,__( ! There is a conjecture in the Physics literature called the scaling hypothesis which suggests 

O I that the large scale asymptotics at the critical temperature should agree with the short 

(Ih ' distance behaviour of the massive scaling theory. A rather detailed version of this scaling 

hypothesis has been proved for the two point function [12], [13]. Modulo the fact that this 

scaling hypothesis can only be literally correct for the scaling limit from above Tc (since 

the critical Ising correlations vanish for an odd number of spins on the lattice and the odd 

r> I scaled correlations from below Tc are non-zero) this hypothesis does suggest a connection 

d ' between the result we consider here and the conjectured result for the critical scaling limit. 

Incidentally, the technique to be explained below should extend to an analysis of the short 

distance asymptotics of the odd correlations in the T <Tc scaling limit. 

In [11] M. Sato, T. Miwa and M. Jimbo showed that the scaled correlations for the 
Ising model have logarithmic derivatives that can be expressed in terms of the solution 
to a monodromy preserving deformation problem. In one version of this connection the 
log derivative is expressed in terms of the solution to non-linear monodromy preserving 
deformation equations. Of more technical interest for us is an alternative connection with 
the linear differential equation which gives rise to the monodromy preserving deformation. 
Much recent work on the asymptotics of non-linear isomonodromic equations exploits in a 
precisely analogous manner the connection with an associated linear problem (see [3] and 
references for an account of modern Riemann-Hilbert techniques in the analysis of such 
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problems) . 

I will begin by recalling the Sato, Miwa, Jimbo characterization [11] for the logarithmic 
derivative of the scaled correlations for the two dimensional Ising model (a more detailed 
proof of this result appears in [8]). The two dimensional Ising model on the lattice T? is 
a classical statistical mechanical system of ±1 spins with ferromagnetic nearest neighbor 
interactions. The Boltzmann weight for a configuration of spins Oi at sites z G A with A a 
finite subset of T? is given by 

where the sum is over nearest neighbors < ij > in A {K > favors nearest neighbor 
alignment and is referred to as the ferromagnetic case) . To be a little more precise if the 
sum is restricted to nearest neighbors < ij > both of which are in A we will refer to open 
boundary conditions. If the nearest neighbors < ij > are allowed one element (say i) on 
the boundary of A (but not in A) and we extend the spin configuration to the boundary 
by setting ai = +1 then we will refer to plus boundary conditions. Write 

< (7,^(7^2 ■■■ai„ >T, (1.1) 

for the thermodynamic (infinite volume) limit of the spin correlation functions at temper- 
ature T > 0. For T less than the critical temperature Tc we consider the limit of the finite 
volume correlations with plus boundary conditions and for T > Tc we consider the limit 
with open boundary conditions (these are the limits discussed in [8]). Next we introduce a 
scaling parameter S > and a temperature T{6) < Tc so that the correlation length (the 
exponential rate of clustering for (1.1)) is S~^ and we define for aj G R^, 

r_(ai,a2,...,a„) = lini , ^ . ^ ' (1-2) 

<5^o < a >^^^^ 

where < a >t=< (Ji >t> is the translation invariant one point function (the sponta- 
neous magnetization). The existence of this limit and the subtleties associated with the 
fact that S~^aj is not always a lattice point in Z^ are discussed in [7]. We will refer to r_ 
as the scaling function for the Ising model from below Tc. We define a scaling function, 
r+, for the Ising model from above Tc by 

/ N V < CT^-iaiCT^-iaa ■ ■ '^(J-ia,! >T((5)* .^„x 

r+(ai,a2,...,an) = hm — — ^ , (1.3) 

where T{6)* is the Kramers- Wannier dual temperature T{6)* > Tc (see [7]). Note that 
in the denominator it is T{S) < Tc which appears and not the dual temperature; this is 
significant since < o" >t= for T > Tc. 

In [11] Sato, Miwa and Jimbo introduced a characterization of (ia log r± (a) which we 
will now describe for t_. For later developments it will be useful to introduce a parameter 
m > which will scale the dependence on the aj in (1.2) and (1.3). For Sato, Miwa and 
Jimbo the functions t± are tau functions for monodromy preserving deformations of the 
Euclidean Dirac equation on R^. The appropriate Dirac operator is. 
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Although we will not be working exclusively with holomorphic functions, the presence of 
d and d in the Dirac operator makes it very convenient to identify R^ with C in the usual 
way, (a;i,a;2) ^ z = xi + 1x2 and of course, 2 = xi — 1x2- Henceforth we will use this 
identification and work over C instead of R^. For brevity we will write f{z) for a function 
of the two real variables xi and X2, even though it is customary to use a notation like 
f{z,z) to avoid the temptation to regard f{z) as a holomorphic function of z. 

Let a = {ai, 02, ... , On} denote a collection of n distinct points in C. Let 7^(C\a) 
denote the 2 fold ramified covering space of C\a with base point ao covering ao G C with 
ao 7^ ttj for j = 1,2, .. .n. Let 7^ denote a simple closed loop based at ao surrounding 
aj but no other point ak for k 7^ j. Let Tj denote the deck transformation of 7^(C\a) 
associated with 7^. 

The differential operator vfiI—0 acts naturally on smooth functions / : 7^(C\a) -^ C^. 
We will call a smooth function w : 7^(C\a) -^ C^ an Ising type solution to the Dirac 
equation if and only if 

(ml — 0) w = 0, 

and w oTj = —w for j = 1, . . .n (that is w changes by a sign "going around" the loop 
7j). An Ising type solution to the Dirac equation has local expansions (essentially Fourier 
expansions) about each of the points tti, given by 

wiz) = Y^ cl{w)wkiz^) + cl*{w)wl{z^), 
kez 



where Zi, = z — a^j, the quantities c^(w) and c^*(w) are complex numbers which we refer to 
as local expansion coefficients, and the functions Wk{z) and w1,{z) are multivalued solutions 
to the Dirac equation defined by (see [11]), 



Wk{z) 



w*k{z) 



Vk{z) 

Vk+l{z) 



Vk+i{z) 
Vk{z) 



■^z'- 



-5r(fc+i) 
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^0{m 



fc+T 



(1.4) 
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2" + 3r(A: + f) 



2'=-5r(fc+i) 
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(1.5) 



Where Vk{z) = e^'^^ "^'^ Ij^_i[mr)^ z = re*^ is the polar form for z. Is is the modified Bessel 
function of order s, and F is the usual gamma function. The functions Wk{z) and wl{z) 



are simultaneously solutions to the Dirac equation (ml — - 



]w 



(*) 



and eigenfunctions 



for the infinitesimal rotation that commutes with the Dirac operator. Observe that the 



functions Wk{z) and wl{z) become less and less locally singular as functions of z as the 
index k increases. Whenever we write down a local expansion that explicitly exhibits only 
a few terms it is understood that the remaining terms + ■ ■ ■ involve less locally singular 
wave functions w^* (i-e., higher values of the index k) than the terms which make an 
explicit appearance. 
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If we define C = 
following conjugation, 
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then it is easy to see that wT is obtained from Wk by the 



wl{z) = Cwk{z). 



The conjugation defined on / : C ^ C^ ^J f ^ f* '■= Cf will have special significance 
for us as it commutes with the Dirac operator and defines a real structure on the space 
of solutions to the Dirac equation. This real structure plays a role in the Ising model 
application we are looking at; observe that because the monodromy multiplier -1 is real 
the conjugation * also defines a real structure on the space of Ising type solutions to the 
Dirac equation. 

Now we recall some results from [11]. Write x ■ y — xiyi + X2y2 for the standard 
bilinear form on C^, so that the standard Hermitian form is given hj x ■ y. If w and v 
are Ising type solutions to the Dirac equation then the inner product u ■ v descends to 
a function on C\a since the monodromy multipliers cancel out. We will say that such a 
solution, w, is in L'^{C) if and only if 

- w ■ w dzdz < oo. 

It is a result of SMJ [11] the space of Ising type solutions to the Dirac equation that are 
in I/^(C) is an n dimensional vector space. The local expansions for such a solution, w, 
are restricted, 

oo 

^(^) = XI Cki'w)'^k{z^) + cl* {w)wl{z,,) for i^ = 1, . . . , n, 

fc=0 

and there is a formula for the inner product of two such L^ solutions in terms of the leading 
local expansion coefficients, 

I f u-vdzdz = -^Y.^^''o*(v) = -^Y.^^''o(v), (1-6) 

Z I (-< Til Til 

^ 1/ 1/ 

which follows from the identities, 

i It 

-u ■ V dzdz = d(u2Vidz) = — d(uiV2dz), (1-7) 

2 mm 

that are a consequence of u and v satisfying the Dirac equation, and the use of Stokes' 
theorem to localize the integrals of the differentials on the right of (1.7) to the singularities 



tti,. Equation (1.6) shows that L^ Ising type solutions are determined uniquely by the low 
order expansion coefficients Cq, z/ = 1, . . .n (or the coefficients Cq* ,1^ = 1, . . . ,n). We will 
now characterize the L^ Ising type solution we are interested in by specifying the behavior 
of its low order expansion coefficients. This solution will satisfy a "reality" condition with 
respect to the conjugation * so it will be convenient to work with a real basis constructed 
from Wk and w^. Define ^Wk = ^{wk + w^) and Qwk = ^{^k — wl). Then for each 
z/ = 1, . . .n there exists a unique Ising type solution Wi, to the Dirac equation which is 
real (W* = Wj^), is in L^, and which has leading order local expansions, 

W^z) = 6^,'^wo{z^) + T^^^woiz^) + ■ ■ ■ , (1.8) 

for z near a^. The coefficient matrix T^^, is real since Wi, is real but as we shall see it is 
not necessary to specify anything further about it to achieve uniqueness (it will turn out 
to be skew symmetric). We will consider the existence and uniqueness for Wi, in a moment 
but first we note that the connection with the scaling function r_ for the Ising model is, 

<i„ log r_ (ma) = ^^cU^A^i.)da^ -c^JW^da^. (1.9) 

V 

This result is due to SMJ [11] (see also the proof of theorem 5.1 in [8]) and will be the 
foundation of the analysis we make of the to — *> limit of r_ [md) that concerns us in this 
paper. Note that in [11] it is zWjy that is taken as the fundamental wave function and this 
alters the formula (1.9) just a little. 

We now recall the construction of the wave function Wj^ and the proof that it is 
uniquely determined by the conditions given above. In [11], SMJ prove that there exists a 
canonical basis {w,^} for the space of 1? Ising type solutions to the Dirac equation. This 
basis is characterized by local expansions with leading order terms, 

^v{z) = SfiiyWoiz^) + C^*{w^)Wq{z^) H , (1.10) 

for z near a^. Applying the conjugation * to (1.10) we get what is called by SMJ the dual 
canonical basis, 

w*(2;) = S^r,WQ{Zf,) + c'f^*{w^)wo{z^) H (1.11) 

The formula for the inner product (1.6) shows that L^ Ising type solutions are uniquely 
determined by the values of the n linear functionals Cq for z/ = 1, . . . , n. Hence we find, 

w.(z) = 5^cr(w.)w:(^), (1.12) 

a 

and 






c^*(w,)c^*(w,) = <5^,. (1.13) 



If we use the formula (1.6) to calculate the inner product of two elments of the canonical 
basis we find, 

i f _ _ 4 * 4 

<w^,w^>=-/ Wf,-w^dzdz = -c(^*(w^) = 5-Co*(w^) (1.14) 
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Now write Cg for the matrix with iju matrix element Co*(w,^) and we can summarize (1.13) 
and (1.14) by saying that Cq is an hermitian symmetric, negative definite matrix with 
CqCq = /. If we now rewrite the focal expansion result for Wi, in terms of wq and Wq 
instead of '^wq and Kwq we find by comparing local expansions that, 






(1.15) 



2 / y\-H-^ ■ --M^^^^^/x- 

Comparing these two results with (1.12) we find, 

c*{T-iI) = {T + iI), (1.16) 

where T denotes the matrix with fii/ matrix element T^j^. We can solve this for T to get 

T = i{cl-I)-\4 + I). (1.17) 

The fact that Cg is negative definite implies that Cg — / is invertible so we can define T by 
(1.17). Since Cg is Hermitian symmetric it follows that T* = —T. Since CgCg = / it follows 
that T — T, that is T is real. Thus T must be (real) skew symmetric as well. With this 
definition for T we can define Wi, by (1.15) and uniqueness follows from the developments 
leading up to (1.17). 

We are interested in computing the m -^ limit of the right hand side of (1.9). In 
order to do this we will first introduce a slightly different characterization of Wi, which 
will introduce the principal player in our understanding of the m ^ limit; this is an 
appropriate Green function for the Dirac operator ml — acting on a suitable space of 
"multivalued" functions (or more precisely it is an explicit formula for the m = limit 
of this Green function). We begin by noting a uniqueness result which will suggest the 
definition of the Green function we are interested in. Suppose that w is an L^ Ising type 
solution to the Dirac equation which has leading order local expansions, 

wiz)=c''{w)^wo{z^) + ---, (1.18) 

for z near a^. Then w is identically 0. To see this observe that for such a w one has 
Cq(w) = Cq*{w) = 2 ■ Substituting this result in the formula for < w,w > obtained 
from (1.6) one finds that < w,w > must be negative unless c^{w) = in which case 
< w,w >= 0. 

Roughly speaking our strategy in dealing with Wi, can now be described as follows. 
We will make a local subtraction (p^, from Wi, which will kill off the term S^,^'<swo in the 
expansion (1.8) for Wu- The difference Wi, — <pi, then satisfies the inhomogenous Dirac 
equation, 

imI-0){W^-(P,) = -{mI-0)<P^:=,U, (1.19) 

and Wjy — (pi, is uniquely characterized by this equation and the fact that it is in I/^(C) 
with local expansions of type (1.18) at each of the points a^. We will demonstrate the 
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convergence of the solution to (1.19) in the hmit to — > with sufficient control to say what 
happens to (1.9) in this limit. 

We will now sketch the existence theory needed to make this possible (we follow the 
existence theory in SMJ [11] quite closely). Let Cq^ {7l{C\a)) denote the space of complex 
valued C°° functions, /, on 7^(C\a) such that f oTj = — /, and such that the projection 
of the support of / in 7^(C\a) onto C\a is compact. Let ?i^° denote the Hilbert space 
completion of C(^(7^(C\a)) with respect to the norm determined by the inner product, 

<f,g >i,o= I j (W{z)dg{z) + m^J{z)g{z)^ dzdz. (1.20) 

Let Ti.^^ denote the Hilbert space completion of Co°(7^(C\a)) with respect to the inner 
product 

<f,g>oX-='^j (W{z)dg{z) + m^J{z)g{z))dzdz. (1.21) 

Now suppose that / G C(f^(7^(C\a)) and consider the linear functional, 

Co~(7^(C\a)) 3g^ '- J W)9i^)dzdz. 

For ?T2 ^ it is clear that this linear functional extends to a continuous linear functional on 
both TC}^ and TiJ^^. Thus by the Riesz representation theorem there exists an F*-^ G 7i^^ 
so that, 

< F''^,g>,,j= '- I J{z)g{z)dzdz, (1.22) 

where i,j = 1, or 0,1. It follows from this last equation that away from the branch points 
F^'^ is a distribution solution to 

(m^ _ A)F''^ = f. 

Elliptic regularity implies that F^'^ is C°° "away from the branch points" (i.e., on the 
covering space 7^(C\a)). Near the branch point z = a^, F^'^ is an L^ branched solution 
to the Helmholtz equation (m^ — A)F'^'^ = and so has a local expansion, 

F^'^iz) = J2 <''''Mz -a.) + Yl b7''Mz - a.), (1-22) 

fc>0 A;>0 

Since dF^'^ is locally in L^ it follows that Oq' ' =0 and since dF^'^ is locally in L^ it 
follows that 6q' ' = for z/ = 1, . . . , n. 

Now we apply this to solve the Dirac equation, 

{m-0)F = f, (1.23) 

7 



where / 



and fj e C^(7^(C\a)). Let F"'i G H'^^ be the solution to (m^-A)F' 



7i' 

/i constructed above and let F^'*^ G Ti^*^ be the similar solution to (m^ — A)F^'° = /2. 
Define, 



(1.24) 



It is easy to check that i^ is a solution to the Dirac equation (1.23) and furthermore that 
because i^*'-' lies in Til^ it follows that the solution F is locally in L^. This is not quite 
the solution to (1.23) which we want. Instead, notice that because F is locally in L^ about 
each point a^, we can subtract a linear combination of the canonical wave functions W;^ 
from F to produce a new solution F to (1.23) with the property that, 

c;^(F) = Ofor z/ = l,...,n (1.25) 

Observe that these boundary conditions uniquely determine an L^ solution to (1.23) since 
the formula (1.6) shows that any L^ solution to the homogenous equation killed by all the 
linear functionals Cq is necessarily 0. Now that we know existence and uniqueness for L^ 
solutions to (1.23) with the boundary conditions (1.25) we can use this enlarge the class of 
suitable boundary conditions to include the conditions relevant to the Ising model. Let A^ 
denote the space of L^ Ising type solutions, w, to (to —0)w — 0. Let J\f denote the image 
of A^ under the map, 

N3w^ {co{w),c*o{w)) G C" © C", 

where Cq{w) denotes the n tuple (cq (tw), . . . , Cq (w)) in C". We know that Af is an 
n dimensional subspace of C^ © C^ and that the intersection of A/" with {0} © C" is {0}. 
Since {0} © C^ is n dimensional it follows that Af and {0} © C" are transverse subspaces 
in C"^ © C"^. Now let X denote any subspace of C" © C^ which is transverse to A/". Then 
we claim that for any / G Cq^ {7l{C\a)) there exists a unique L^ Ising type solution, F, 
to {m — 0)F = /, with, 

(co(F),c^(F))GJ. (1.26) 

Existence is simple. We know that the differential equation (to — 0)F = f has an L^ Ising 
type solution with co{F) = 0. By transversality there is a unique element w & N so that, 

0(Bc*{F)-{co{w),c*{w))eI. 

It is clear that F — w is a, solution satisfying the boundary condition (1.26). Uniqueness is 
obvious. 

Now we introduce the boundary condition that is relevant to the Ising model. Hence- 
forth, let X denote the subspace of C" © C" which consists of all vectors of the form {v, v) 
with V G C". This corresponds to the boundary condition, 

co(F) = cS(F), (1.27) 

which the reader should note is the "complexification" of condition (1.18). Clearly X is n 
dimensional and the formula (1.6) for the inner product of L^ solutions shows that any 



element w E M which has boundary values in X has an L^ norm which is negative unless 
the linear functional cq vanishes on it. This is enough to show that the function must be 
and hence that X fl A/" = {0}. Thus X is transverse to M and there is an existence and 
uniqueness result for solutions to (1.23) satisfying (1.27). This is the complexification of 
the result we really want. If / is a real function, that is, 



f = r 



fl 



then it is easy to see that the solution, F, satisfying (1.27) must also be real (i.e., F* = F). 
In this case both co{F) and Cq{F) will be real numbers. 

It is straightforward to describe the subtraction (p^, that we make from Wj/. Let fu{x) 
denote a real valued C°° function on C which is identically 1 in a small ball centered at 
tti, and which vanishes outside a slightly larger ball about a^, but still small enough so that 
ipiy is zero in a neighborhood of all the other points a^ for /i 7^ z/. Define 

Consulting (1.8) we see that Wi, — (f)i, will have local expansions at each point a^ of type 
(1.18). Write SW^, := m2(Wi, — (piy). Then using the fact that W^, and ^wq both satisfy 
the massive Dirac equation we find. 



{m-0)SW^ 



-2av7,. 

-2dip^ 



m^Qwoizi,) := fu 



(1.28) 



The scale factor ma is introduced here so that following limit exits, 

1 



lim m^^wo{z,y) - 

m^O \/2'Ki 



(1.29) 



In terms of 5Wi, the coefficient md{ {Wi, 
by, 



mdi{W.) 



which appears in the formula for d log r_ is given 




- (<5W.)i (ee^ 



-^de^. 



(1.30) 



The function [5^^^]^ is the first component of SWu and the integral on the right of (1.30) 
calculates a (half integer) Fourier coefficient of this function on the circle of radius e about 
a^j. Since 



lim — 

m^O /i 



1 



(me) 



it will suffice for our purposes to control the m ^ convergence of dWi, in LP{C^{a,y)) 
for any p > 1 and all u. Here C^{a,y) is the circle of radius e about tti,. A device that will 
be useful for this purpose will be to "excise" small neighborhoods of the branch points 
tti,. Let Di, denote a small open disk centered at tti, and suppose that the radii of these 
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disks are small enough so that they are mutually disjoint. We also want D^, chosen small 
enough so that the ip^, is identically 1 on D,^. The function f^, will then have its support 
outside the union of the D^y. Write C^, — dD^ for the clockwise oriented circle that is the 
boundary of D^j. Let D^o denote a disk centered at ao with a radius that is large enough 
so that all the disks Di, are contained in the interior of -Doo- We will write W^ for the 
complement in C of Doo • Write Coo for the counterclockwise oriented circle dDoo ■ Let V 
denote the bounded domain which is the complement of the union of the closed disks Dj^ 
for z/ = 1 , . . . n in Doo • We now rephrase the problem we need to solve to find 5 Wjy . We 
wish to find a solution to 

ioY X eV with the restriction of 5Wv to C^ for ^u = 1, . . . , n equal to the boundary value 
of a solution to the Dirac equation in D^ with local expansion of type (1.18) (we will refer 

to this subspace of boundary values as W/i"^ ) . We also require that the restriction of SWu 
to Coo belongs to the space of boundary values of solutions to the Dirac equation which 
are in L^(D^) (we refer to this subspace as W^ ). Of course, we really need to work in 
a two fold covering space for T) and two fold covers of the circles Cq, but for simplicity of 
exposition in the remainder of this note we will ignore this (the details of the argument we 
only sketch here will appear in a paper that is still in preparation). Our principal tool in 
understanding the to ^ limit of the solution to this problem is a "guess" for the m = 
Green function. Here is the guess, 



1 

Go{x,y) -- 



Y.u'^^nAx)vy{y) 9{x,y ) 

9{x,y) ^„u^{x)v^{y) 



Airi 

where 

u^{x) = {x-a^} 2 I I tf_^ 



9{x,y) = > c(e)^ 

, , y — X 

with e = (ei, e2, . . . , tn) and each e,^ = ±4. Also 



a 



v ) 






E 



e,., 



c(6) := 
and 



]\u<u\^^^-(^^ 



|2ef.ei. 



|2e^ey 



{x-a^Y^ 



Vj,{x) = (x - ttj.) 2 ^ c(e) JJ^ 

|e(i/)|=0,e, = i i^^i^ 



{at, - a^)^M 



A possible surprise in this guess are the "chiral symmetry breaking terms" on the 
diagonal in Gq. The m = Dirac operator is purely off diagonal and one might expect 
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the same for its Green function. However, the boundary conditions we are interested in 
force the off diagonal terms even in the m = hmit. Indeed the principal ingredient that 
went into this guess was setting the boundary conditions for the m = limit equal to the 
limits, W}, , of the subspaces W},"^' as to — *> 0. It is not difficult to check that Gq is a 
Green function for —0 on T) essentially using the fact that -^- is a Green function for d. 

X y 

As a first approximation to dWi, we invert m — according to the following scheme, 

(m - 0)-^ = -0-^{l - m$-^)-^ = Go(l + mGo)"^ (1.31) 

It is a classical result [5] that the Cauchy kernel -^- defines a bounded operator on L'p{T>) 
for all p > 2. This result carries over without difficulty to Go{x,y) (modulo the difference 
that we should be working on a covering of V). The Neuman series for (1 + mGo)~^ then 
converges in the strong operator topology on L'p{T>) and we can use this to show that 
(1.31) does invert m — for all small enough m. Our ffist approximation to SWi, we take 
to be 

-^iW, :=Go(l + mGo)"V.- 

This solves the appropriate inhomogeneous Dirac equation but the boundary values are 
in the subspaces Wj) instead of Wu . The to — > limit of 5iWi, is straightforward to 
compute. 

The Calderon projector [2] associated with the Green function Go^x^y) maps the 

direct sum (BW^ onto ©Wj) along the subspace, N, of boundary values of solutions 
to 01/) = 0. For small to one can show that this projection determines an isomorphism of 
(BWu with ©Wj) with a complementary projection on N which is 0{m). The fact that 
this is true for the Green function Go is precisely what makes our choice of Go suitable 
in this application. Thus we can adjust diWi, by an element 62^1, G N which is 0(171) so 
that SiW,y + 52W1, has boundary values in Wo' . This is not quite the solution we want 
but since 52'Wv G N it follows that 

(m - 0) {5iW, + 52W,) = /, + TO(52W,. 

Thus the adjustment we need to get SWu from 5iWv + 52'Wv is given by the solution to 
the Dirac equation with right hand side 17162)^^1, = 0{m^) and boundary values in (BW},^ . 
We don't have very detailed control of the solution to this problem but because the right 
hand side is O(to^) a simple apriori estimate shows that this last adjustment is 0{m) as 
well. Making use of the explicit form for the m — > limit of 61W1, we are now in a position 
to calculate the limiting form of the logarithmic derivative of r_. We find, 

lim dlogT- = -dlog V Ff I^a* " ^1^1^'^'"- (1-32) 

\e\=OfJ.<i' 

This is the Luther-Peschel formula for the 2n point function. 
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